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1. Introduction 

At the Large Hadron Collider (LHC) as well as at the Tevatron, electroweak di-boson 
production in association with a hard jet represents an important class of processes, of 
either signal or background character in various searches for Standard Model (SM) and 
beyond. The rates are large, especially at the LHC by accessing the gluon density at small 
momentum fraction, and next-to-leading order (NLO) QCD corrections have turned out 
sizable in a series of recent publications [1-5] providing WW-l-jet, ZZ+jet, and W^7-|-jet 
production at NLO QCD precision at hadron colliders. In this paper, we supplement NLO 
QCD precision to 

pp, pp ^ 3 leptons -|- + jet -|- X , 

i.e. to W^Z-|-jet production including full leptonic decays. We give cross sections for 
LHC and Tevatron collisions and also discuss the corrections' phase-space dependence by 
investigating differential correction factors at the LHC. 

To verify our results, special care is devoted to independent numerically stable imple- 
mentations of the processes, yielding two independent fully-flexible Monte Carlo programs, 
based on different approaches. 

We organize this work in the following way: In section ^ we provide an overview of the 
two programs we have employed for the numerical results of this paper, to which section ^ 
is devoted. We first focus on on-shell production at the Tevatron and the LHC in sections 



3.2 and p^, respectively, and then move on to discuss the differential impact of the QCD 



corrections for the LHC setup in more phenomenological detail in section 3^. Section 
concludes with a summary of the work presented in this paper. 
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Figure 1: Representative Feynman graph contributing to pp,pp — >■ Pge^/i^/i^ + X. The crosses 
mark points where the 7, Z ^ fJ-^ fJ-" decay topology can be inserted. Indicated is also the exchange 
of a virtual gluon, which gives rise to self-energy, triangle, box, and pentagon topologies. Not 
shown are topologies that result from closed fermion loops, non-abelian graphs, and real emission 
topologies. All other subprocesses can be recovered from the shown graph by flavour summation 
and/or crossing. 

2. Details of the calculation 

We invoke a dedicated system of checks and balances to validate our calculation. In par- 
ticular, we have compared two different Monte Carlo implementations, relying on distinct 
approaches. The comparisons involve cross checks at the amplitude level for a fixed phase- 
space point as well as comparisons of integrated cross sections. In both cases we find 
agreement on the level of double-precision accuracy or agreement within the statistical 
errors, respectively, for different choices of renormalization and factorization scales and 
cuts. 

Our independent calculations are based on the approaches of [2-4], where the used 
methods have already been described in detail. Hence, we limit ourselves to the bare 
necessities to make this paper self-consistent and refer the interested reader to the above 
publications for more details. 

Program 1 We generalize the NLO QCD calculation of W7-|-jet production (including 
leptonic decays) [4] to WZ-|-jet production. The leading order (LO) matrix elements at 
0{a^as), cf. figure H] are calculated with Helas routines [6] generated with MadGraph 
[7]. Our phase-space implementation is based on routines readily present in the Vbfnlo 
suite [8], which were already applied in the context of NLO QCD vector-boson- fusion 
WZ-|-2jets production, including leptonic decays, in various scenarios [9,10]. Nonetheless 
integrated results for the different subprocesses have been checked against Sherpa [11], 
yielding agreement within statistical errors on per-mill level. 

The virtual corrections are combined to groups that include all loop diagrams derived 
from a born-level configuration, i.e. all self-energy, triangle, box and pentagon corrections 
to a quark line with three attached gauge bosons, which are computed as effective decay 
currents in case of the electroweak bosons, are combined to a single numerical routine. 
This leaves a set of universal building blocks, which were already appropriately assembled 
to determine the one-loop contribution to W7-|-jet production (cf. [4] for details on the 
verification of the implementation against an independent approach). The generalization 
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to WZ+jet becomes trivial by replacing the photon polarization vector by the effective Z 
decay current multiplying the appropriate coupling. These building blocks are set up us- 
ing in-house routines within the framework of FeynC ALC [12] and FeynArts [13] . They 
invoke the Passarino-Veltman reduction [14] up to boxes and the Denner-Dittmaier reduc- 
tion [15] for pentagons. The remaining fermionic loop corrections are derived via algebraic 
calculations using FeynCalc, which is subsequently processed to Fortran routines with 
in-house routines. The scalar integrals that arc not already present in the Vbfnlo frame- 
work are supplemented from the Ellis-Zander ighi library [16]. All effective decay currents 
are evaluated by means of Helas routines generated with MadGraph, which are modified 
to fit our purpose of calculating the one-loop amplitude. 

To speed up the numerical implementation of the numerous subprocesses that show 
up as part of the real emission, we computed the real emission matrix element using the 
spinor hclicity formalism of [17]. We store intermediate numerical results common to all 
subprocesses and re-use them whenever possible. All matrix elements have been checked 
explicitly against code generated from MadGraph. 

The infrared (IR) singularities are subtracted applying the dipole subtraction of [18], 
while the corresponding LO matrix elements and currents for the subtraction kinematics 
are computed using Helas routines. We also apply necessary bookkeeping in order not 
to waste computing time. The IR poles of the virtual amplitude are cancelled against the 
real-emission ones algebraically, and we perform the integration of finite coUinear terms as 
part of the real-emission integration by appropriately mapping the born-type configuration 
as done in [19]. 

The code will become publicly available with an upcoming update of Vbfnlo. 

Program 2 We proceed essentially in the same way as in the calculation of WW-l-jet, 
which is discussed in some detail in Refs. [2,3]. All LO helicity amplitudes are calculated 
by application of the Weyl-van-der-Waerden formalism (as described in Ref. [20]). In this 

approach, the implementation of the gauge-boson decays can be easily realized by replacing 
the polarization bispinors of the gauge bosons with expressions containing the currents of 
the decay leptons. Allowing for off-shell gauge bosons and still respecting gauge invariance 
requires the inclusion of diagrams which do not contain two simultaneously resonant gauge- 
boson propagators. However, all diagrams of this kind can be constructed from W-|-jet- 
production amplitudes by replacing the W-boson polarization bispinor with an appropriate 
expression describing its decay into 4 leptons. Additionally, diagrams with intermediate 
photons instead of Z bosons have to be taken into account. It is worth noting that the 
described replacements are exactly the same for the LO and for all components of the NLO 
QCD calculation — in other words, the bispinor replacements are universal. In particular, 
no new types of loop diagrams show up in the virtual corrections. 

Again, the dipole subtraction formalism of [18] is applied to rearrange the IR diver- 
gences between real and virtual corrections at NLO QCD. 

The loop diagrams and amplitudes are generated by FeynArts 3.4 [13] and then 
further manipulated with FormCalc 6.0 [21] to automatically produce Fortran code. 
The whole reduction of tensor to scalar integrals is done with the help of the LoopTools 
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library [21], which also employs the Denner-Dittmaier method [15] for the 5-point tensor 
integrals, Passarino-Veltman [14] reduction for the lower-point tensors, and the FF package 
[22, 23] for the evaluation of regular scalar integrals. The IR (soft and collinear) singular 
3- and 4-point integrals in dimensional regularization are linked to this library as in the 
WW+jet calculation. Again, the explicit results of Ref. [24] for the vertex and of Ref. [25] 
for the box integrals (with appropriate analytical continuations) are taken. Actually the 
FormCalc package assumes a four-dimensional regularization scheme for IR divergences, 
i.e. rational terms of IR origin are neglected by FormCalc. However, in Ref. [26] it was 
generally shown that such rational terms consistently cancel if UV and IR divergences 
are properly separated. Thus we could use the algebraic result of FormCalc for the 
unrenormalized amplitudes without any modification, apart from supplementing the needed 
IR-singular scalar integrals. 

To receive the real-correction matrix elements we also employ the Weyl-van-der- 
Waerden formalism. The dipoles needed to cancel the divergences in the respective sub- 
processes are automatically generated from the born-level helicity amplitudes. To achieve 
numerical stability on a high-accuracy level, the phase-space integration is performed by a 
multi-channel Monte Carlo integrator [27] with weight optimization [28], which has been 
written in CH — h and checked in detail in the calculation of WW-|-jct. Additional channels 
basing on dipolc kinematics arc automatically included to improve the efficiency of the 
integration of the dipole-subtracted real-emission matrix elements. When the full calcula- 
tion with off-shell gauge bosons is considered, only channels according to doubly-resonant 
diagrams are included, which turns out to provide already sufficient numerical stability. 

Both numerical programs account for finite width effects of the electroweak gauge 
bosons (when considering their leptonic decays) with a fixed-width scheme, which is also 
the scheme used by MadGraph: while we calculate with Brcit-Wigner propagators of the 
W and Z bosons we keep the weak mixing angle real. To justify this approach, which 
breaks gauge invariance, we compared the results to a calculation performed applying the 
complex-mass scheme [29] in one of the programs. We find an agreement on the per-mill 
level between the two calculations, so the effect of gauge-invariance breaking turns out to 
be sufficiently small to be ignored here. 

3. Numerical results 

Throughout, we use CTEQ6M parton distributions [30] at NLO, and the CTEQ6L1 set at 
LO. We choose Mz = 91.1876 GcV, Mw = 80.425 GeV, and Gp = 1.16637 x 10"^ GeV"^ 
as electroweak input parameters and derive the electromagnetic coupling a and the weak 
mixing angle sin^^ via SM tree-level relations. The LO and NLO running of as are 
determined by a^^(Mz) = 0.130 and a^^'^{Mz) = 0.118 for five active flavors, respectively. 

The center-of-mass energy is fixed to ^/s = 14 TcV for LHC and ^/s = 1.96 TeV 
for Tevatron collisions, respectively. We consider both on-shell production of the elec- 
troweak bosons and their decays to distinct species of light leptons, e. g. W~ — )■ e^Pg and 
Z — )• jjL^ 1 treating these leptons as massless. The CKM matrix is taken to be diagonal, 



-4- 



and we neglect bottom contributions throughout because they are numerically negligible 
anyway and can even be further suppressed by b-tagging. To be more precise, we neglect 
the — finite and negligibly small — contribution from real correction minus subtraction terms 
if external bottom quarks are involved. In the fermionic quark loops and, correspondingly, 
the /-operator, we keep all six quark flavours. A non-diagonal CKM matrix decreases our 
LHC results only at the per-mill level because gluon-induced subprocesses dominate the 
cross section. In case of a Cabibbo-like block-diagonal CKM matrix, the contribution from 
this subset of subprocesses is not affected if all light quarks are summed over. The cor- 
rection for the Tevatron cross section is about 3% due to the dominance of quark-induced 
subprocesses. These corrections are well below the residual scale dependence at NLO QCD. 
The final-state partons are recombined to massless jets via the algorithm of [31] with res- 
olution parameter R = 1.0. Other jet algorithms, like the kT algorithm of Ref. [32] have 
also been implemented in program 1. 

3.1 Event selection 

To analyze the impact of the NLO QCD corrections on the total production cross sections 
at both the Tevatron and the LHC, we apply a rather inclusive set of cuts. In case of 
on-shell W and Z boson production, the jets are required to have a transverse momentum 
of 



PTjet > 50 GeV , (3.1) 
which is the only selection criterion we impose for the calculation of the cross sections of 



Sees. 3.2 and In case of included leptonic decays of the electroweak bosons, we account 



for finite jet-detection coverage by requiring the jets to have rapidities 

|7?,| <4.5 



(3.2) 



in addition to the cut on pxjet of Eq. (3J.). However, with the used cut value for PTjet, 
the additional effect of this cut is completely negligible at the given collider energies. All 
leptons are required to lie in 



\Vi\ < 2.5 



(3.3) 



with transverse momenta of 



PT/ > 25 GeV . (3.4) 
The overall missing transverse momentum is chosen to be 

> 25 GeV . (3.5) 
The leptons have to be separated in the azimuthal-angle-pseudorapidity plane by 

Ree> = (A(/.|, + Ar/|^,)^^^ > 0.2 , (3.6) 
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and for their separation from observable jets we choose 

Rij>OA. (3.7) 

It is customary from a theoretical point of view to study also the cross sections' behaviour 
with a veto applied to the second resolved jet in the context of mono-jet production: it was 
shown [1-5] that this additional veto yields highly stabilized NLO exclusive cross sections 
in the context of WW+jet, ZZ+jet, and W7+jet production. Indeed, given the similarities 
of the processes from a QCD point of view, identical properties for WZ+jet production are 
evident. We will discuss the phenomenological problems of the additional jet veto in detail 



in section 3.4. 



3.2 Production cross sections at the Tevatron 

For on-shell W^Z+jet production from Tevatron collisions with the additional requirement 
of ( |3.lD , we compute a total inclusive and exclusive WZ+jet on-shell cross sections and K 
factors 

NLO 

K = -^ (3.8) 

of 

f^Sci°(W^Z + jet) = (139.01 ± 0.10) fb [K = 1.209] , (3.9) 
f^^xc?(W^Z + jet) = (129.40 ± 0.10) fb [K = 1.125] , (3.10) 

which are dominated by qQ induced processes due to the relatively large momentum frac- 
tion of the incoming partons x ~ 0.2 at LO, which we infer from the Monte Carlo simulation. 
The total correction of 21% with respect to LO at the central scale is sizable. Nonetheless, 
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including the leptonic decays decreases the cross sections to phenomenologically subdom- 
inant size, unless the transverse momentum requirement for jets is reduced substantially. 
We therefore limit ourselves to quoting total on-shell production rates at the Tevatron and 
focus on differential distributions at the LHC only. 

A lower bound on the scale uncertainties of the cross sections can be inferred e. g. 
from varying the fixed renormalization and factorization scales by a factor of two around 
the central value //r = /ip = Afwj cf. figure Doing so, the LO approximation exhibits a 
scale variation of 31% which is decreased to 9% by including NLO-inclusive precision. 

3.3 Production cross sections at the LHC 

Turning to the more energetic LHC collisions, we find a completely different situation 
compared to the Tevatron. The proton is typically probed at much lower momentum 
fractions x ~ 0.02 at LO (as the Monte Carlo calculation shows), so that the g^-induced 
initial states dominate the total rate. The total NLO-inclusive cross sections are 

f^Sci°(W"Z + jet) = (7.495 ± 0.008) pb [K(W-Z + jet) = 1.298] , (3.11) 
f^Sci°(W+Z + jet) = (12.061 ± 0.013) pb [K(W+Z + jet) = 1.260] , (3.12) 

for a scale choice /iR = = Mw- 

In close analogy to [1-5], the di-jet contribution re-introduces a substantial dependence 
on the renormalization scale ^r via the dominating (/^-induced channels. This becomes 
apparent by checking the variation for several scales intrinsic to the total cross section in 
figure ^. 

The cross sections' qualitative scaling behaviour does not depend on the choice of the 
intrinsic scale, and the characteristic increase of the NLO-inclusive cross sections at small 
scales /i = /iR = /ip reflects the renormalization scale dependence of the di-jet contribution, 
which is a leading order-Og contribution to our NLO computation. 

This /Ur dependence of di-jet contributions can be effectively buffered by imposing an 
additional veto on events with two resolved jets, which gives rise to total NLO-exclusive 
rates for ;^r = /^f = M^j of 

o-^xc?(W~Z + jet) = (4.981 ± 0.009) pb [i^(W-Z + jet) = 0.862] , (3.13) 
f^^xc?(W+Z + jet) = (7.831 ± 0.014) pb [E:(W+Z + jet) = 0.818] . (3.14) 

Varying again /iR = /ip by a factor two around the central values in figure |3| amounts 
to scale uncertainties of 8% (W~Z+jet) and 8% (W+Z+jet) of the total inclusive cross 
sections at the LHC. For the vetoed sample, the scale dependence is reduced to about 5% 
for W~Z+jet, and 6% for W^Z+jet production. 

The improved perturbative stability of the exclusive cross sections should be inter- 
preted with caution. While jet-vetoing is a straightforward exercise in the context of fixed- 
order Monte Carlo calculations, its phenomenological consequences are generally highly 
delicate, both from the theoretical and the experimental side. The small total correction 
along with the stability against variations of /iR = /ip of the exclusive cross sections should 
therefore not be misinterpreted as a guideline to stable LHC predictions per se, but as a 
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Figure 3: Scale variation for on-sliell W^Z+jet (upper plots) and W^Z+jet (lower plots) produc- 
tion at the LHC. The identified renormalization and factorization scales /xr, — fip — ^ are varied 
with respect to the fixed scale Mw, the maximum px.jet, and the invariant mass towz- The plots 
on the left show cross sections at LO, those on the right at NLO QCD accuracy. 



significant perturbative improvement of WZ+jet production up to the specified threshold 
value of PT,jet- This is visible in the differential jet-pT distribution of Fig. where the 
uncertainty band is particularly narrow for small transverse momenta. Whether this ad- 
ditional jet veto gives rise to a sufficiently stable theoretical approximation in the sense 
of an experimentally applicable strategy, does highly dependent on the phenomenological 
question we ask, i.e. the phase-space region we are interested in. Additional jet radiation, 
as can already be inferred from figure ^, is kinematically unsuppressed to large extent, es- 
pecially when considering hard events with large transverse momenta. Vetoing additional 
radiation in a region of phase-space where it becomes likely is crucial to the flat scale 
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Figure 4: LO, NLO-inclusive and NLO-exclusive differential distributions of the transverse mo- 
mentum of the hardest jet and its pseudorapidity for W+Z+jct production including leptonic decays. 
The differential if-factor band corresponds to varying /iR = /ip by a factor two around the central 
scale in the NLO distribution only. 



dependence of the exclusive cross section. We discuss this in more detail in Sec. 3.4. 
3.4 Differential distributions at the LHC 

We now turn to the effect of QCD corrections to the full processes pp ^ 3 leptons + 
+ jet + X. Including the leptonic decays with the selection criteria quoted in section 3A 
yields the cross sections given in table |l| at //r = //p = -^w for the inclusive and the vetoed 
sample, respectively. In table || we additionally give a precision comparison of the cross 
sections calculated with our two programs described in Sec. ^. 

The differences of W~Z-|-jet compared to W^Z+jet production are predominantly due 
to the different parton distribution functions of the incoming partons in the dominating 
subprocesses. In particular, there are no initial-state up-quarks involved in W~Z-|-jet at 
LO, but in W+Z+jet. 



Program 1 



Program 2 



a 



NLO 

incl, decay 



NLO 



cr. 



cr. 



incl, decay 
NLO 

excl, decay 
NLO 

excl, decay 



(W-Z+jet) 7.4592 [48] fb 7.4628 [63] fb 
(W+Z+jet) 11.129[10]fb 11.1286 [47] fb 



(W-Z+jet) 4.6721 [62] fb 
(W+Z+jet) 6.6900 [92] fb 



4.6663 [64] fb 
6.6816 [49] fb 



Table 1: Comparison of the numerical results from both of our programs to verify their excellent 
statistical agreement on the per-mill level for pp 3 leptons + + jet + X. 
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Figure 5: LO, NLO-inclusive and NLO-exclusive differential distributions of the maximum and 
minimum lepton transverse momentum for W~Z+jet production. The differential iiT- factor band 
corresponds to varying /ir = yUp by a factor two around the central scale in the NLO distribution 
only. 



From figures ^-|^ we uncover a substantial observable-specific phase-space dependence 
of the QCD corrections. While additional jet radiation gives sizable contributions to the 
maximum-pTjet distribution at large values, the LO approximation considerably overesti- 
mates the NLO-exclusive findings. Additionally, for inclusive events, the jets tend to be 
more central due to the extra hard jet emission, which occurs central at small rapidity 
differences with respect to the other reconstructed jet. 

The harder inclusive jets balance against a softer inclusive lepton-px spectrum, figure |5|, 
and additional QCD radiation. Apart from this significant deviation of the leptonic dis- 
tributions from the total cT'^^^/cx'^^-rescaling, the bulk of the leptonic observables receive 
only minor differential distortions compared to rescaled-LO when including inclusive-NLO 
precision. Representatively, we show the maximum and minimum charged- lepton pseudo- 
rapidity in figure ^ and the tri-lepton invariant mass 

"^feptons = (Pe- + Pf^+ + P^" )^ (^-l^) 

and the transverse WZ cluster mass in figure ^. The transverse cluster mass 

"iT.ciuster = Y^m2 (e-/i+/x- ) +p|(e-/i+/x-) + I^tI) - (PT(e"/^+^") +#t)^ ' (3-16) 

is a convenient observable to observe production of additional charged heavy bosons [10,33] 
from a beyond-the-SM sector via Jacobian peaks. 

The exclusive distributions, even though improved perturbative stability is suggested 
from the decreased scale dependence of the total cross sections, which can be observed in 
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Figure 6: LO, NLO-inclusive and NLO-exclusive differential distributions of the maximum and 
minimum lepton's pseudorapidity for W^Z+jet production. The differential if-factor band cor- 
responds to varying /iR = /ip by a factor two around the central scale in the NLO distribution 
only. 



figure ^, exhibit large uncertainties, especially in the tails of the pT distributions. Here 
additional jet radiation becomes probable as can be seen from the maximum jet-px dis- 
tribution in figure ^. The improved NLO stability of the exclusive sample shows up as 
perturbative improvement almost exclusively around the threshold region. For the phase- 
space regions characterized by larger values of pT, applying the additional jet veto does 
not yield a stable result anymore — at least in the chosen setup. As shown in figure |5|, 
perturbative control over the exclusive production cross section is already lost at scales 
of about 100 GeV while the inclusive differential cross section turns out to be reasonably 
stable. Obviously the jet-veto with a fixed pT threshold, although hinting at appealing 
properties by the exclusive cross sections' flat scale variations, does not easily give rise to 
a more reliable cross section prediction within the given order of perturbation theory. 



4. Summary 

We have computed NLO QCD cross sections and differential distributions for W^Z pro- 
duction in association with a hadronic jet at hadron colliders. The calculation has been 
extended to full leptonic final states at the LHC, where they are well-observable. We find 
the total QCD corrections to be sizeable at both the Tevatron and the LHC. At the same 
time they show strong phase-space dependencies in hadronic, semi-hadronic, and especially 
in transverse momentum distributions. Hence, QCD modifications should be taken into 
account in every phenomenological study that employs these processes. 
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Figure 7: LO, NLO-inclusive and NLO-exclusive differential distributions of the tri-lepton in- 
variant mass TOicptons and the transverse cluster mass for W^Z+jet production. The differential 
if-factor band corresponds to varying /ir = /ip by a factor two around the central scale in the NLO 
distribution only. 

In addition, we demonstrate that the superficial perturbative improvement for exclusive 
production at the LHC, which is also observed in the various other massive di-boson+jet 
production cross sections [1~5], does not give rise to perturbatively stable predictions once 
the additional-jet veto's impact on the large pT region is taken into account. The exclusive 
production's reduced scale variation therefore expresses NLO stability for a part of our 
calculation which actually is given to NLO precision. This discrimination between NLO 
QCD one-jet and LO di-jet contributions, which is inherent to fixed-order calculations, can 
not be carried over to experimental strategies in a straightforward way, and, hence, does 
not easily give rise to a phenomenologically applicable method. 
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